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ABSTRACT 


Low  frequency  active  sonar  perfonnance  in  shallow  water  is  trflen  limiied  by 
reverberation.  Reverbenutioo  modeling  in  shallow  water  has  been  difficult  due  to  the 
conqdexity  of  die  mult^di  acoustic  piq^gation  problem  inherent  in  shallow 
environments,  hi  August  1992,  a  shallow  water,  low-firequency  reveibecatitm 
measumnent  was  made  in  the  Barents  Sea  utilizing  explosive  "signal,  underwater  sound" 
(SUS)  charges  as  sound  sources  and  a  Iti^lonent  vertical  hydrophcme  amy  as  the 
receiver.  The  objectives  of  diis  diesis  were  to  analyze  die  leverberaticsi  data  from  this 
experiment,  compare  several  theories  which  haw  been  proposed  to  model  reverberation, 
and  (tetermine  the  reverberant  characteristics  (rf  die  regicxi.  The  three-dimensional 
Hamiltonian  Acoustic  Ray-tracing  Program  for  die  Ocean  (HARPO)  was  used  as  the 
primary  jm^gation  modeling  ttxd.  The  tetxqxxal  signal  processing  ccmsisted  of  a  short- 
time  Fourier  transform  spectral  estimation  method  applied  to  data  from  a  single 
hydrc^hcHie.  Chiqmian’s  source  spectrum  model  was  used.  Reverberation  models  based 
rni  Lambert’s  law  and  (mmidirectional  backscattering  theory  were  ccmipared.  Lamben’s 


law  was  found  to  be  mcxe  applicable  in  the  Barents  Sea.  A  statistical  analysis  was 
performed  on  broadband  and  narrowband  hydrc^riione  data  showing  that  reverberation  in 
the  Barents  Sea  possesses  Gaussian  properties. 
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L  INTRODUCTION 


A.  MOTIVATION 

As  the  worid  pdiitical  climate  oontinoes  to  shift,  the  new  naval  battlefield  has  been 
extended  into  littonl  regions  (OlCeefe,  1992).  Prosecution  of  a  new  breed  of  coastal 
diesel  subimrines  has  become  a  top  priority  for  ensuring  U.S.  farces  maintain  battlespace 
dominance  in  diese  regioas.  The  develoimient  oi  low  frequency  active  acoustic  (LFAA) 
systrais  by  several  nations  worldwide,  including  the  United  States,  may  provute  new, 
more  effective  metiiods  by  which  such  forces  can  be  combatted.  Such  systems  represent 
a  double-edged  sword,  however,  because  tiieir  potential  effectiveness  in  the  detection  and 
prosecution  of  diesel  submarines  operating  in  littoral  regions  implies  a  threat  to  the 
tactical  security  of  U.S.  submarines.  Hence,  performance  predictions  for  LFAA  systems 
have  become  a  priority  for  the  Navy. 

The  performance  of  LFAA  systems  in  shallow  water  is  often  limited  by 
reverberation  from  the  bottom.  Reverberation  modeling  is  a  coiiq)lex  task,  for  it  involves 
neatly  every  aspect  of  underwater  acoustics  research.  It  incocpOTates  environmental 
mcmitoring,  prcqiagation  modeling,  underwater  explosion  theory,  anay  analysis,  and 
scattering  theory  into  a  single  field  of  research.  Because  the  problem  is  so  complex, 
reverberation  modeling  has  been  historically  difficult  due  to  the  lack  of  adequate 
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processilif  power.  Hence,  even  diou^  die  need  for  eccurate  predictions  has  exined  for 
many  years,  only  recently  has  it  been  possible  m  adequately  model  the  problem. 

In  the  past,  much  research  efifon  has  been  expended  in  the  study  of  reverberatimi 
in  deep  water.  Most  of  the  experiments  have  involved  relatively  simple  prc^agadon  and 
scattering  models  (McCammon,  1993).  They  have  mainly  used  point  sources  and  single 
hydrophone  receivers  to  analyze  enei^gy  which  traveled  cm  a  single  padi  to  a  scatterer  and 
then  returned  on  that  same  path.  This  lessened  the  computational  load  significantly,  but 
limited  the  range  fw  which  the  analysis  was  accurate.  Nonetheless,  in  deep  water,  diese 
relatively  sinqile  models  were  adequate  to  produce  relatively  accurate  performance 
predicticms  of  sonar  systems. 

In  shallow  water,  due  to  the  multipath  acoustic  propagation  taking  place,  such 
models  were  not  adequate.  Hence,  a  more  accurate,  and  necessarily  more  complicated 
modeling  technique  was  needed.  This  need  motivated  the  inclusion  of  a  low  frequency 
reverberation  study  as  part  of  the  Barents  Sea  Polar  Front  Experiment  (BSPFEX).  In  this 
experiment,  a  vertical  array  and  explosive  "signal,  underwater  sound"  (SUS)  charges  were 
utilized  to  ctxiduct  low  frequency  reverberation  measurements.  These  measurements  then 
provided  baseline  data  which  have  been  used  in  this  thesis  to  test  the  validity  of 
reverberatitm  prediction  models. 

The  experiment  also  provided  a  valuable  qiportuniQ'  to  investigate  the  statistical 
pn^ierties  of  reverberation  in  the  Barents  Sea.  In  a  shallow  water  environment  with  a 
rigid  bottom  such  as  the  BSPFEX  experiment  area,  the  highly  multipath  nature  of  acoustic 
propagation  creates  a  reverberant  signal  which  is  a  combination  of  many  random 
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processes,  each  of  which  may  or  may  not  be  Gaussian.  By  the  central  limit  theorem, 
however,  tl»  sum  of  these  signals  anriving  at  a  single  hydrophone  should  possess  a 
nmrnal  Gaussian  distribution.  This  hypothesis  was  examined  as  part  d  this  BSPFEX  data 
analysis  for  both  narrowband  and  broadband  «gnal«- 

B.  OBJECTIVES  AND  APPROACH 

The  main  dmists  of  this  thesis  can  be  summed  up  in  three  objectives.  The  first  was 
to  examine  two  commonly  used  laws  for  reverberation  prediction  as  applied  to  the 
Barents  Sea,  namely  Lambert’s  law  and  omnidirectional  scattering.  The  second  objective 
was  to  estimate  the  environmentally-dependent  reverberation  nxxlel  parameter  |i  for  the 
Barents  Sea  and  to  determine  its  frequency  (tependence.  The  third  objective  was  to 
perform  a  statistical  analysis  of  the  BSPFEX  data  to  determine  if  the  characteristics  of 
reverberation  in  the  region  did  in  fact  possess  Gaussian  properties. 

Considering  the  oceanogiiqrhic  data  available  fiom  unclassified  sources  and  the 
computing  resources  on  hand,  the  fiuUowing  methods/inodels  were  used  for  theoretical 
predictions  and  data  analysis  in  this  thesis.  Narrowband  data  analysis  was  performed 
using  shOTt-time  Fourier  transforms  for  all  power  spectral  density  estimates.  Statistical 
analysis  was  performed  on  broadband  and  narrowband  data  using  normalized  data 
segments  and  histograms.  Propagation  modeling  was  perfOTmed  utilizing  ray  theory  as 
implemented  in  the  Hamiltonian  Acoustic  Ray-tracing  Program  for  the  Ocean  (HARPO), 
with  post-processing  performed  on  the  HARPO  ouq)ot  to  determine  the  effects  of 
spreading  and  boundary  interactions.  Geometric  effects  were  accounted  for  using  a  ray- 


3 


tate  spVMdbf  alforitlan.  Boimdary  inienctioiis  were  accounted  for  using  the  Rayleigh 


reflection  coefficient  and  Eckait  rough  surface  scattering  function.  The  prt^gation  due 
a  point  source  was  qtfnoximated  using  four  ray  launch  azimuths  ("4  x  2D”  estimate). 
Backscattering  functions  developed  by  Lambert,  Tolstoy  and  Clay,  and  McCanunon  were 
investigated.  A  SUS  source  model  devel(q>ed  by  Quqiman  (1988)  was  used.  Finally,  the 
propagatitMi  model,  backscattering  function,  and  source  model  were  combined  to  produce 
a  narrowband  reverboaticm  envelope  prediction  for  the  Barents  Sea.  By  comparing 
theoretical  incdictions  for  the  Barents  Sea  with  the  experimental  measurements,  the 
validity  of  each  reverberation  theory  was  tested. 

The  most  significant  choice  of  the  die  models/methods  listed  above  was  that  of  the 
propagation  model,  for  it  had  the  greatest  effect  on  the  level  of  the  predicted  signal  and 
it  determined  which  models  or  methods  were  i^ipropriate  for  boundary  interactions  and 
scattering.  All  routines  developed  for  this  diesis  were  programmed  in  the  Matlab* 
environment  and  implemented  on  an  HP  9000  series  735  computer  (TAC  3),  a  general 
purpose  UNDC  based  machine. 


C.  OVERVIEW  OF  THESIS 

The  approach  described  above  is  summarized  in  Figure  1.1.  Environmental  data 
were  gathered  fnxn  unclassified  sources,  including  bathymetric  and  geologic  data  from 
recent  publications  and  sound  speed  field  data  determined  during  the  BSPFEX.  HARPO 
was  used  to  trace  rays  on  several  azimuths  from  the  array  location  in  order  to  estimate 
the  three-dimensional  propagation  effects  which  determine  the  transmission  paths  between 
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were  con^Mred  to  determine  the  validity  of  the  scattering  functions  used.  The  statistical 
analysis  of  the  BSPFEX  data  was  performed  as  previously  described  and  compared  with 
a  ntvmal  Guissian  distribution  to  verify  its  statistical  properties  agreed. 

In  Cluq>ter  II,  details  of  the  reverberaticm  experiment  in  the  BSPFEX  are  presented. 
The  discussion  includes  a  brief  discussion  of  oceanographic  and  geologic  features  of  the 
Barents  Sea,  a  description  of  the  SUS  charges  used  as  sound  sources,  and  the  construction 
and  operation  of  the  vertical  array.  This  is  followed  in  Qiapter  HI  with  a  discussion  of 
reverberation  thetxy,  both  in  general  terms  and  specific  theories  which  have  been 
developed.  Theories  to  be  presented  include  Lambert’s  Law,  Tolstoy  and  Clay’s  low 
frequency  scattering  function,  and  McCamnxm’s  low  frequency,  low  grazing  angle 
backscattering  strength  model.  Chiqrter  IV  presents  a  method  for  predicting  bottom 
reverberation  using  HARPO,  a  source  model,  and  a  backscattering  strength  nnodel.  It  also 
presents  the  results  of  the  BSPFEX  and  conqjares  them  to  the  predicted  levels.  It  ends 
with  an  analysis  of  the  statistical  properties  of  the  reverberation  measured  during  the 
experiment.  Hnally,  Chapter  V  summarizes  the  work  completed  and  discusses  lessons 
learned.  Advice  regarding  future  work  to  refine  die  present  study  will  be  provided. 
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our  potm  fouroe  (SUS),  botioin  locatioiis,  and  a  hydrophone.  A  SUS  source  spectrum 
model  from  recent  rnearch  was  used  as  a  second  input  to  the  reverberation  prediction. 
The  last  factor  needed  was  the  bottom  scattering  strength  function  Sk(f»6t»0).  Scattering 
functions  suggested  by  Lambert,  Tolstoy  and  Clay,  and  McCammon  were  considered. 
With  these  inputs,  tt»  reverberation  level  present  at  the  array  was  predicted. 


Figure  1.1:  Reverberation  level  measurement  and  prediction  flowchart. 


The  measurement  portion  of  the  problem  was  addressed  in  a  somewhat  more 
straightforward  manner.  A  spectral  analysis  subroutine  was  applied  to  the  hydrophone 
data  to  produce  a  time  series  of  reverberation  level  at  specific  frequencies.  Frequencies 
of  SO,  100,  200,  and  400  Hz  were  considered  in  this  thesis.  These  finequencies  have  no 
special  significance  (i.e.  to  an  existing  LFAA  system),  other  than  the  fact  that  they  cover 
the  array’s  50  to  500  Hz  bandwidth  in  octave  steps.  The  measured  and  predicted  levels 
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n.  BARENTS  SEA  EXPERIMENT 


A.  EXPERIMENT  OBJECTIVES 

The  Barents  Sea  Polar  Front  E]q)eriinent  (BSPFEX)  was  ccmducted  in  August  of 
1992.  The  eiqwriment  was  a  jdnt  effort  hetwem  the  Naval  Postgraduate  School  (NFS). 
Woods  Hole  Oceant^ny^  Institute  (WHOD,  and  the  Science  Applications  International 
Corporation  (SAIQ.  The  principal  instigators  for  the  experiment  were  Professors 
Ching-Sang  Chiu,  James  R  Miller,  and  Robert  Bourke  from  NPS,  Dr.  James  F.  Lynch 
frxMn  WHOL  and  Dr.  Robin  Muench  frtmi  SAIC  The  principal  engineers  for  the 
develqmient  and  deployment  of  the  vertical  hydit^hmie  array  system  were  Mr.  Keith  von 
der  Heydt  and  Mr.  John  Kenqi  bam  WHOL  The  objectives  of  this  experiment  as 
outlined  by  the  Barents  Sea  Polar  Front  Groiq)  (1992)  were: 

1.  Provide  a  detailed  physical  descriptitm  of  the  polar  front 

2.  Enhance  the  understanding  of  dynamics  of  the  front  including  frontogenesis 
and  its  influence  on  regional  o^anogn^hic  processes. 

3.  Assess  the  ability  of  acoustic  tcnnognqthy  to  define  frontal  and  associated 
mesoscale  features. 

4.  Provide  improved  documentation  of  shallow  water  acoustic  propagation  in 
this  region  and  the  effect  of  the  environment  on  acoustic  Anti-Submarine 
Warfare  (ASW)  operations. 


7 


LONGITUDE. 

Figure  2.1:  The  Barents  Sea. 

The  reverberation  experinaent  was  designed  to  meet  part  of  the  fourth  goal  by  studying 
low  frequency  sound  propagation  and  scattaring  in  shallow  water.  (McLaughlin,  1993) 

B.  EXI^RIMENT  LOCATION 

The  experiment  was  performed  in  tire  Barents  Sea  100  km  east  of  Bear  Island,  as 
shown  in  Hgure  2.1.  The  Barents  Sea  ftxms  an  epicontinental  sea  (1.3  million  km^), 
bounded  by  the  Arctic  Ocean  to  the  north,  the  Svalbard  archipelago  and  the  Norwegian- 


Batfiymtry  In  Expartmant  Aiaa  (|n  ntalara) 


Longituda  Eaat 

Figure  22:  Bathymetry  in  the  BSPFEX  area. 

Greenland  sea  to  die  west,  the  Fennoscandian  shield  in  the  south  and  Franz  Josef  Land 
and  Novaya  Zendya  to  die  east  It  is  characterized  by  northeast-southwest  trending  basins 
300  to  500  meters  in  depth,  with  an  average  depth  of  230  meters  over  the  entire  region. 
It  is  thus  signiBcandy  deeper  than  most  of  the  present  day  high  Arctic  shallow  shelves 
outsitte  North  America,  Nathem  Europe  and  Northem  Asia  (which  are  generally  10  to 
60  meters  deep).  The  greater  depth  is  most  lilmly  a  resptmse  to  repeated  glaciations  in 
the  Late  Cenozoic,  leaving  only  a  thin  sediment  cover  above  the  Mesozoic  and  Paleozoic 
bedrock  (typically  5  to  .15  meters  thick)  (Norsk  Polarinstitutt,  1987). 
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The  majoricy  of  the  Barents  Sea  is  covered  with  a  thin  sediment  layer,  though  large 
sediment  distiibutkms  are  present  in  water  depths  exceeding  300  meters  in  the  western 
part  oi  die  mi^  troughs  Bj^rn4yrenna  and  StotQordrenna  and  exceed  SOO  meters  in 
thickness  near  the  shelf  edge.  In  the  experiment  area,  the  sediment  thickness  varies  from 
nearly  0  msec  (exposed  bedrock)  in  die  northern  portion  to  approximately  SO  msec  near 
the  southern  extent  of  the  area.  The  majority  of  the  sediment  is  composed  of  stiff  pebbly 
mud,  dll  and/or  glacimnarine  deposits  overrun  and  reworked  by  glacial  activity  (Norsk 
Polarinstitutt,  1987). 


Sound  SpMd  PreO*  Nmt  Sm  Vwtieal  Um  Airay 


Figure  23:  Sound  speed  inofile  near  the  vertical  line  array. 

Near  the  experiment  area,  the  oceanography  is  dominated  by  the  Barents  Sea  Polar 
FitMiL  It  is  a  stable  feature  which  is  caused  by  circulation  in  the  region  (Emblidge,  1991) 
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and  which  passes  directly  through  the  experiment  area  ap|»t)ximately  S  km  north  of  the 
receiving  array.  The  {Aysical  oceanognq;^y  oi  the  region  plays  a  large  role  in  controlling 
sound  propagation.  The  bottom,  as  previously  described,  is  rough  in  some  areas  and 
ransists  of  a  thin  sediment  overiying  hard  rock.  The  badiymetry  (Figure  2.2)  and  sound 
speed  are  highly  range  dependent  The  shallow  water  and  generally  downward  refiracdng 
sound  velocity  imrfile  cause  many  surfsce  and  bottom  interacticms  in  the  acoustic 
pit^Migation. 

The  sound  speed  field  was  observed  to  vary  significantly  with  range.  The  sound 
speed  profile  measured  at  die  vertical  line  array  is  shown  in  Figure  2.3  and  is  a  typical 
Barents  Sea  profile  south  of  the  front  for  the  summer.  The  sound  speed  gradient  is 
basically  negative  from  surfrce  to  bottom  widi  a  shallow  mixed  layer  region.  The 
negative  gradient  tends  to  bend  rays  toward  the  bottom,  enhancing  tiie  region’s  already 
strong  botKxn  reverberation  characteristics  and  minimizing  surface  reverberant  effects. 

C.  EQUIPMENT  DESCRIPnON 

The  locations  for  the  vertical  line  array  and  the  SUS  charges  which  were  processed 
are  shown  in  Figure  2.4  and  are  listed  in  Table  2.1.  A  description  of  each  follows. 

1.  SUS  Charges 

A  total  of  24  SUS  charges  were  expended  during  the  experiment  Eighteen 
of  the  charges  were  delivered  by  a  P-3  Orion  aircraft  which  overflew  the  data  recording 
ship  on  two  courses,  242*T  and  351*T,  dropping  9  SUS  charges  on  each  run.  The  SUS 
charges  dropped  by  the  P-3  were  U.S.  Navy  Mark  64-0  exercise  communication  SUS, 
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Figwre  14:  Device  locations  for  leverbeiation  ^qperimenL 
whidi  conaan  a  1.1  lb  ^.50  kg)  TNT  main  diazge  and  were  set  to  eiqdode  at  60  ft  (18 
meters).  The  charges  have  a  nominal  broadband  source  level  of  263  dB  re  l|iPa  @  1 
meter  (U.S.  Navy,  1989).  All  of  the  air-dropped  charges  e;q)loded  and  were  recorded,  but 
the  60  foot  dettmadon  dqxh  placed  the  charges  well  widiin  the  mixed  layer.  Although 
they  were  easily  heard  by  the  array,  the  reverberation  created  by  the  charges  was  miniinal 
and  die  high  propagation  loss  between  the  detonadon  points  and  the  array  severely 
attenuated  die  signal  Hence,  the  recordinp  of  die  air-dropped  SUS  are  of  little  use  for 
this  particular  research. 

The  leoiainigg  six  charges  were  dirown  by  handcrff  the  back  of  die  lecoiding 
ship.  One  of  the  six  was  a  dud,  aixi  one  was  not  recorded.  The  ship-dropped  SUS 
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TABLE  2.1:  DEVICE  LOCATIONS 


Device 

latinide 

Longitude 

Distance  to 

Amy 

Vertical  Line 
Amy 

74*  19.1994’N 

23*  32.2957’E 

N/A 

SUS  3 

74*  19.06*N 

23*34.70*B 

1232  meters 

SUS  4 

74*  18J4’N 

23*  33.68'E 

1405  meters 

SUS  5 

74*  18.97*N 

23*  34.72’E 

1287  meters 

SUS  6 

74*  18J9’N 

23*  34.20’E 

1479  meters 

charges  used  were  U.S.  Navy  Maik  82  special  purpose  SUS,  which  contain  a  1.8  lb  (0.82 
kg)  TNT  main  charge  and  are  set  to  explode  at  a  depth  of  300  ft  (91.4  meters).  The 
charges  have  a  nominal  broadband  source  level  of  274  dB  re  IpPa  @  1  meter  (U.S.  Navy, 
1989).  They  were  dropped  at  an  average  range  of  1350  meters  from  the  array  and 
produced  a  revetberadon  signal  which  is  distinguishable  above  the  background  noise  level 
for  approximately  SO  sectxids. 

Qualitatively  speaking,  Mk  82  SUS  provide  an  inq;>ulsive  signal  which, 
although  non-linear  in  the  immediate  vicinity  of  the  source,  serves  as  a  very  close 
estimate  of  a  point  source  with  an  iiq)ulse-like  pressure  signature.  A  time  series  of  one 
such  expk>sion  is  shown  in  Hgure  2.5.  As  noted  in  the  following  paragraphs,  the  signal 
would  be  even  nxne  in^ulsive  in  qipeaiance  if  it  had  not  been  for  non-linearities  in  the 
receiving  array.  A  more  complete  discussicm  of  underwater  explosion  theory  is  discussed 
in  Chapter  m.  Reverberation  Theory. 
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Figure  2.6:  Venical  line  anay  and  RF  Ethemet 
1992) 


buoy.  (V(m  der  Heydt.  et  al.. 


In  the  tomognqiriiy  nqperiment,  the  sign  (+/-)  of  the  signal  made  little 
difference,  since  die  (xily  significant  data  was  the  travel  time  of  the  signals.  In  the 
reverberation  experiment,  it  was  found  tiutt  tiie  received  signals  were  consistently  180° 
out  of  phase  with  dieoiy  due  to  the  details  of  the  data  acquisitirai  system.  This  shows  up 
most  vividly  in  the  initial  blast  waveform,  which  should  be  initially  an  overpressure  in 
tile  ground  wave,  then  tiie  waterborne  shock  wave,  again  positive-going,  followed  by 
several  positive-going  bubble  pulses.  After  the  initial  arrivals,  arrivals  from  other 
multipaths,  including  surfsce  reflections  (which  do  undergo  a  180°  phase  shift  and  hence 
would  be  initially  negative-going)  come  in  but  are  corticated  by  the  bubble  pulses  of 
tiie  previous  arrivals.  The  corrected  signal  is  shown  in  Figure  2.7.  In  all  data  analysis. 


IS 


8U8  3  Ftwiw  18  Nydropheiw  Vofeao* 


-0^  -0.1 
TifM(sac) 

Figure  2.7:  Initial  received  waveform  firom  Mk  82  SUS. 
the  polarity  of  the  received  rignal  will  be  reversed  in  order  to  account  for  this  feature. 

An  additional  undocumented  feature  of  die  vertical  line  array  was  discovered 


during  analysis  of  the  SUS  explosion  time  series.  Either  the  hydn^hones  on  the  array 
or  the  analog-ttHiiigital  (A/D)  converter  in  the  buoy  showed  a  non-linear  response  at 
levels  exceeding  ~160  dB  (1  volt),  witii  the  signal  entirely  clipped  above  165  dB  (1.8 
vcrits).  Hgure  2.8  illustrates  this  clearly  during  die  first  half  second  of  the  received  blast 
The  signal  stould  have  been  a  sharp  increase  to  a  peak  (at  ~19S  dB),  followed  by  a  near 
linear  decrease  in  amplitude  (in  decibels).  This  would  correspond  to  an  exponential 
decrease  in  pressure  with  time.  Instead,  the  signal  remains  at  die  maximum  phone  output 
level  until  the  initial  shock  wave  and  first  bubble  pulses  pass.  The  A/D  converter  had 
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Mk  82  SUS  Uv«l  at  1 Ion.,  Rtotiwd  Phon*  0  of  Array 


Tima  aflar  axplosion  (sac) 


Figure  2J:  Reveibentioa  kvd  pkx  showing  clipping  of  hydn^hone  ouq)ut 

been  qiedficaUy  designed  to  show  a  linear  response  over  -5  volts  to  S  volts  hydrophone 
ouq>uL  Therefore,  die  frixmes  themselves  were  likely  the  source  of  the  non-linear 
behavior.  To  be  conservative,  it  is  assumed  that  the  phones  display  a  ntm-linear 
anqtiitude  response  over  ISO  dB  re  1  uPa,  as  is  illustrated  gnqihically  in  Figure  2.9. 

The  in^lication  of  this  on  the  reverberation  experiment  is  that  approximately 
the  first  second  of  the  SUS  explosion  data  is  unusable.  This  is  unfortunate,  since  any 
time-domain  methods  of  reverberation  require  a  time  series  as  the  input  signal  to  the 
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lewteiiloB  xpiarioii.  Bamiw  of  dds.  ^>ec«Ba  levdt  of  SUS  duBges  fiom  recent 


reaench  have  been  used  in  all  data  analysis  and  prediction  routines. 


Mppws  nyorapnonv  owsviQf  wvw 


Figure  2J9:  Apparent  hydiO|dione  lesptxise  curve. 


m.  REVERBERATION  THEORY 


A.  GENERAL  REVERBERATION  THEORY 

Any  discontinuity  in  the  [rfiysical  pr(^)erties  of  a  medium  tend  to  inteicept  and 
reradiate  a  portion  of  the  acoustic  energy  incident  upon  them  (Urick,  1983).  This 
reradiadon  oi  energy  is  called  scatterings  and  the  sum  total  of  scattering  from  all 
scattcners  in  the  medium  to  a  receiver  is  called  reverberation.  These  discontinuities  can 
be  due  to  any  type  of  inhomogeneity  present  in  die  sea.  Particles  of  inhomogeneous 
material  (dust,  fish,  bubbles,  etc...),  the  sea  surface,  the  bottmn,  and  changes  in  sound 
velocity  in  the  medium  away  from  a  perfect  isovelodty  structure  (vertical  or  horiztmtal 
effects)  all  cmitribute  as  scatterers. 

Reverberaticm  can  be  generally  categorized  into  two  broad  areas:  Volume  and 
surface.  Volume  reverberation  is  due  to  energy  which  has  undergone  scattering  with 
some  inhmngeneity  which  is  suspended  in  die  water  volume  and  then  has  propagated  back 
to  the  desired  receiver  location.  Exanqiles  of  scatterers  which  contribute  to  volume 
reverberaticm  include  marine  life,  suspended  sediment  particles,  and  bubbles.  Surface 
reverberation  is  due  to  energy  which  has  undergone  scattering  with  some  inhomogeneity 
in  the  form  of  a  flat  plane.  The  two  most  commonly  considered  types  are  sea-surface 
reverberaticm  and  bonom  reverberaticm.  Several  ^pes  of  volume  reverberation,  however, 
can  be  effectively  modeled  with  an  equivalent  surface  reverberant  layer.  Perhaps  the  most 
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wett-known  exanqple  of  this  method  is  the  scxalled  deep  scattering  layer,  which  is  in  fact 
a  thick  layer  of  marine  life  which  tends  to  vary  in  depth  depending  on  the  time  of  day 
and  die  seastm.  Due  to  its  plane-like  structure,  it  is  typically  more  efficient  and  equally 
accurate  to  model  its  effect  using  a  surface  of  equivalent  scattering  strength  located  in  the 
proximity  of  the  layer  volume. 

Due  to  the  statistical  nature  of  the  large  number  of  scatterers  contributing  to  the 
total  reverberation  level  at  any  one  time,  reverberation  is  usually  discussed  in  decibels 
referenced  to  1  pPascal  (dB  re  IpPa)  instead  of  describing  it  as  a  time  series  in  terms  of 
pressure  for  a  particular  medium.  As  such,  the  most  general  form  of  the  bistatic 
reverberation  equation  is: 


RL  •  SL  ~  "  ’^KOtttnr-netivtr  * 

SL  is  the  source  level  of  the  projects  used  for  the  active  sonar  system  and  is  typically 
given  in  decibels  referenced  to  1  pPascal  at  one  meter  (dB  re  IpPa  @  1  m).  For  sources 
which  are  beamfonned,  this  term  may  be  different  depending  on  the  angle  from  which 
the  energy  leaves  the  source.  SL  may  be  given  as  die  total  level  over  a  band  of 
frequencies  or  as  a  narrowband  power  spectral  density  referenced  to  a  ^  Hertz  bandwidth 
(units  of  dB  re  1  pPaVHz  @  1  m).  The  two  TL  terms  represent  the  transmission 
(propagation)  loss  from  the  source  to  the  scattering  inhonxigeneity  (scatterer)  and  from 
the  scatterer  to  the  intended  receiver  of  the  sonar.  Energy  traveling  from  the  source  to 
the  scatterer  on  one  "path"  will,  in  most  cases,  be  scattered  back  to  the  receiver  along 
many  paths  connecting  the  scatterer  and  the  receiver.  For  this  reason,  the  propagation 
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modeling  which  is  used  both  in  prediction  and  data  analysis  methods  can  be  the  most 
significant  choice  that  is  made  in  the  reverberation  nxxleling  process. 

The  last  term,  (reverberant  target  strength)  is  the  ratio  (usually  given  in 
decibels)  of  iiwident  power  arriving  at  the  scatterer  on  one  path  to  the  power  which  is 
retransmitted  along  a  return  path  to  the  receiver.  This  term  is  dependent  tm  both  the 
charactoistics  of  the  volume  or  surface  causing  the  scattering  and  the  total  area  or  volume 
of  the  scatterer.  There  are  two  possible  ways  in  which  such  interactions  may  be 
predicted.  The  first,  discrete  backscatter  prediction  (Medwin,  1981,  Clay  et  al.,  1977,  and 
Dyer  et  al.,  1993),  involves  a  detailed  knowledge  of  the  scattering  surface.  Each  scatterer 
("wedge",  "plane  strip",  or  "facet",  respectively)  must  be  identified  and  its  effect  on 
incoming  energy  quantized  in  order  to  effectively  predict  the  resulting  backscattered  field. 
Every  interaction  of  incoming  energy  in  a  given  prt^gation  mechanism  (eigenray  or 
mode)  produces  an  impulse  in  the  backscattered  field  which  is  both  weighted  by  a 
magnitude  and  phase  term  and  shifted  to  a  time  which  is  the  sum  of  the  travel  times  from 
the  source  to  the  scatterer,  and  back  to  the  receiver.  This  method  is  analagous  to 
estimating  the  impulse  response  of  the  ocean  as  if  it  were  a  linear,  time-invariant  filter 
(Ziomek,1993).  The  total  output  reverberation  signal  is  thus  determined  by  a  convolution 
of  the  input  signal  (in  the  time  domain)  with  the  impulse  response  (also  in  the  time 
domain).  As  indicated  by  the  discussion,  these  methods  require  such  precise  knowledge 
of  the  ocean  envinnunent,  which  is  in  fact  time-varying,  as  to  prohibit  accurate  prediction. 
The  distinct  advantage  of  this  iq>proach,  howevo*,  is  that  it  produces  a  signal  estimate 
which  is  corrected  for  phase  differences  due  to  propagation  and  scattering,  and  hence  may 
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idnw  finer  (tail  (liigher  imotoioa)  than  do  the  "envelope"  methods  discussed  below. 
Mediods  which  empk^  discrHe  backscatter  piedictioa  have  been  developed  by  Medwin 
(using  Biot  and  Tolstoy’s  "wedge"  theoiy,  1981),  Qay,  et  al.  (using  "plane  strips",  1977), 
and  Dyer,  et  al.  (using  "fiicets",  1993). 

The  seccmd  method,  leveiberation  envelope  prediction,  ennuis  estimating  die 
scattering  strength  (ISy  or  SJ  of  the  surface  causing  the  leverfaeraticm  depending  on  the 
angle  or  wave  number  of  die  incoming  and  (wtgoing  energy  and  the  type  of  surface 
and/oe  its  estimated  roughness.  A  great  deal  of  research  effort  has  been  expended  to 
determine  scattering  strength  empirically  fat  every  area  in  the  world  over  a  wide  range 
of  fiequencies.  Sy  and  are  given  in  decibels  per  unit  volume  or  area  (Sy  is  given  in 
decibels  referenced  to  1  metei^  (dB  re  m^.  in  decibels  referenced  to  1  meter’  (dB  re 
m’) ).  They  are  in  turn  muldplied  by  the  volume  or  area  of  the  scatterer  (by  adding  10 
/og(Volume)  or  10  log(Area)  in  decibels)  to  give  the  reverberant  target  strength  TSg  of 
the  entile  scatterer.  Only  the  anqilitude  level  (i.e.  not  the  phase)  of  the  reverberant  signal 
is  predicted,  hence  the  name  envelope  prediction.  The  equations  which  generate  the 
reverberant  target  strength  are  thus: 

75it  »  5|r  10  logiVobtme)  2) 

-  54  10  log^rea) 

In  the  BSPFEX,  the  feequencies  of  concern  are  determined  by  the  bandwidth  of  our 
receiving  array,  SO  to  SOO  Hz.  In  this  freqi^ncy  range  the  effects  of  volume  reverberation 
can  be  assumed  negligible  after  1  second  following  the  source  impulse  (NRL  8721, 1991). 
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For  reasons  discussed  at  length  in  Chapter  n.  the  first  second  of  data  will  not  be 
analyzed.  Hence,  volume  reverberation  is  not  a  concern. 

This  thesis  will  be  mainly  ctmoeroed  with  surface  reverberation  caused  by  the  sea 
b(»t(»n,  since  the  majority  of  the  reverberation  in  the  BSPFEX  is  assumed  to  be  due  to 
bottom  backscatter.  This  assuii4)tion  is  based  on  the  ooeanogn^hy  in  tire  Barents  Sea  and 
the  frequencies  which  are  under  investigation  and  the  Barents  Sea  environnrent  The 
effect  due  to  the  oceanography  stems  fixxn  the  fact  that  the  sound  speed  profiles  in  the 
Barents  Sea  are  generally  downward  refracting  with  a  strong  mixed  layer,  both  of  which 
minimize  sea  surface  reverberatitxi  and  enhance  bottom  reverberation. 

In  general,  for  a  surface  to  be  effective  at  reradiating  a  particular  frequency  of 
sound,  it  must  have  features  whose  size  are  (xt  the  order  of  a  wavelength  of  sound  at  that 
frequency.  Ftnr  example,  a  sea  surface  with  tms  waveheight  of  0.5  meter  will  be  an 
effective  scatterer  for  10  KHz  sound,  but  a  weak  scatteier  at  best  for  100  Hz  sound. 
Since  the  seas  were  relatively  calm  during  the  experiment  Qess  than  0.1  meter  rms  wave 
height),  surface  backscatter  can  be  assumed  negligible  compared  to  the  effea  of  the 
bottom.  Experimental  data  nreasured  by  Chapman  and  Harris  suppon  this  assumption  at 
a  frequency  of  500  Hz,  placing  the  surface  backscattering  strength  an  order  of  magnitude 
lower  than  our  bottom  backscattering  strength,  with  the  trend  serving  to  minimize  the 
effect  of  surface  reverberation  further  at  lower  frequencies  (Urick,  1983). 

An  additional  effect  of  the  botunn  is  subbottom  scattering,  the  teradiation  or 
reflection  of  eirergy  fiom  inh(Mix)geneities  below  the  water/sedinrent  interface.  This  effect 
can  also  be  described  as  bottom  volume  scattering.  Subbottom  reverberation  is  not  a 
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m  dke  Baieflti  Sea,  liowever,  due  to  die  extremely  thin  eedtinntt  layer  described 


in  C3i^piBr  II  and  the  "fast”  rock  bottom  immediaiely  below  the  sediment 
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Figure  3.1:  SUS  3  ground  wave  and  wateibome  arrivals. 
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To  verify  this,  the  speed  of  sound  in  the  bottom  can  be  calculated  using  the  SUS 
data.  At  the  leading  edge  of  the  reoei^^  signal  is  the  ground  wove,  eneigy  which  has 
entered  the  bottenn  (at  the  critical  angle  for  the  least  time  path),  propagated  through  the 
bottom  to  the  array,  and  then  reentered  the  water  column  (again  at  the  critical  angle)  to 
be  received  at  the  array.  The  ground  wave  and  waterborne  arrivals  for  one  of  the  SUS 
charges  is  shown  in  Figure  3.1.  The  process  by  which  the  difference  in  travel  times  is 


used  is  diq>layed  pqikicaUy  in  Figure  3J2.  The  speed  of  sound  was  found  to  be  greater 
than  3000  tn^  in  the  bosom,  verifying  that  the  subbonom  scattering  is  insignificant 
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Figure  3,2:  Use  of  ground  wave  to  detennine  sound  speed  in  the  bottom. 

An  additkxial  effea  is  shown  in  Hgure  3.1  immediately  before  the  waterborne 
energy  airives.  The  small  positive-going  spike  ^[iproximately  20  msec,  prior  to  the  direct 
path  blast  arrival  is  most  likely  energy  wbidi  was  trapped  in  the  diin  sediment  layer. 
Using  dw  travel  time  difiference  between  it  and  the  waterborne  energy,  the  sound  speed 
of  the  sediment  layer  is  calculated  as  *  1668  m/s.  This  number  can  dien  be  used 
along  with  the  finquency  content  of  the  trapped  energy  .to  detennine  die  sediment 
thickness.  The  domiruuit  frequency  present  in  the  energy  tr^ped  within  the  sediment 
"waveguide"  is  286  Hz  (accounting  for  die  energy  sdll  propagating  in  die  rock  basement 
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coatfaWBat  with  udNilaied  dati  fbr  die  legioii.  The  signihcance  oi  the  sediment  sound 
speed  and  thickness  lay  in  the  fact  diat  the  mergy  trapped  within  the  sediment  layer  is 
minimal,  and  the  sound  $peed  of  die  sediment  is  such  dut  the  layer  is  almost  acoustically 
tranqpaitnL 

Hence,  widi  sea-surfime  leveiberatioo  and  volume  reverberation  (bodi  in  die  water 
column  and  in  the  bottom)  insignificant  wUh  respect  to  bottom  reverberation,  a  discussitm 
of  bottom  backscattering  theciy  is  warranted. 

B.  BOTTOM  SCATTERING  THEORY 

The  bottom  is  an  effiscdve  reflector  and  scatterer  of  sound  and  acts  to  redistribute 
a  portion  of  die  sound  incident  upon  it  back  into  die  ocean  above  it  (Urick,  1983).  This 
redistribution  of  energy  is  illustrated  in  Figure  3.3,  which  shows  the  effecdve  "beam 
pattern"  of  die  sea  bottom  for  reradiating  incident  sound  energy.  Fmir  cases  are  shown, 
illustrating  the  four  general  cases  which  occur  in  bottom  scattering.  The  first  and  second 
general  cases  (cases  (a)  and  (b))  describe  bottoms  which  have  a  high  impedance  contrast 
with  die  water  volume,  where  die  acoustic  impedance  of  a  medium  is  defined  as  the 
product  of  die  medium’s  density  p  and  sound  speed  c.  Any  sufficiendy  hard  bottom  tits 
diis  definidon,  be  it  rocky  or  hard  sand.  Note  that  almost  all  of  the  incident  energy  is 
redistributed  in  the  water  volume  above,  with  roughness  determining  the  relative  strength 
of  die  specular  reflection.  Cases  (c)  and  (d)  show  the  scattered  field  for  bottoms  which 
have  low  impedance  contrast  with  the  water  column.  Bottoms  of  this  type  are  sometimes 
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High  Zb/Z«v.  Smooth  bottom  High  Zh/Zw,  Rough  bottom 


(c)  (d) 


Figure  3  J:  Directicxial  patterns  of  the  scattered  sound  firom  the  bottom  for 
different  conditions  of  toughness  and  inq)edance  (Za^pc)  omtrast.  (Urick,  1983) 

called  acousticalfy  transparent  in  that  most  energy  is  prc^gated  into  the  botunn,  leaving 

relatively  little  to  be  distributed  back  into  the  water  volume  above.  Examples  of  these 

bottoms  are  soft  mud  or  silQr  bottoms  with  a  thick  sediment  layer.  Note  that  bottom 

rou^mess  is  still  a  dominant  factor  in  the  magnitude  in  bodi  the  specular  and  scattered 

directions. 
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This  brings  iq»  an  iMerewim  Act  Althongb  die  type  of  bottom  sediment  deposits 
may  serve  as  a  first-ctn  means  of  classifying  bottoms  in  tenns  of  acoustic  backscattering, 
the  roughness  of  the  sea  bottom  appem  to  be  the  dominant  factor  in  deiennining  die 
backscattering  characteristics.  Just  as  wind  speed  has  for  many  years  been  used  as  an 
indicator  of  surfoce  loughnen  and  bidible  drasity  near  the  surface,  so  is  the  bottom  type 
merely  an  indicator  of  bottom  rooghness.  (Urick,  1983) 

Many  attempts  to  quantue  and  madiematically  model  bottom  backscattering  are 
documented.  Of  dieae.  three  theories  of  bottom  scattering  will  be  considered  (as 
discussed  by  McCammon,  1993): 

1.  Composite  Roughness. 

2.  Lambert’s  Law. 

3.  Omnidirecticmal  Scattering. 

L  Composite  Rou^mess 

The  first  dieoiy  is  that  of  composite  roughness,  which  McCammem  describes 
as,  "the  most  rigorous  theoretical  ^iproach"  (McCammon,  1993),  since  it  can  be  derived 
themetically  using  the  Helmholtz  integral  equation.  Coiiqiosite  toughness  theory  is 
composed  of  two  fypes  of  scattering  mechanisms:  facet  scatter  and  Bragg  scatter. 

a.  Facet  Scatter 

For  backscatter  due  to  incident  angles  gteater  than  70**  grazing,  and 
forward  scatter  due  to  all  incident  angles,  facet  scatter  theory  is  applicable  (McCammon, 
1993,  and  Tobtoy,  et  aL,  1987).  It  assumes  that  the  bottom  is  composed  of  "facets"  or 
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small  plane-like  smicnires  that  only  ocularly  scatter  ct^ierent  and  inccrfwrent  sound 
energy.  Hdice,  in  order  for  a  fleet  to  scatter  energy  from  a  given  raypath  into  a  scattered 
path,  the  facet  must  be  tilled  sufficiently  to  provide  a  q)ecular  leflecdon.  These  types  of 
models  require  information  regarding  die  bottom’s  rms  rou^mess  height  and  the  spatial 
correlation  length  of  duu  rou^mess.  Most  assume  the  surface's  roughness  is  Gaussian. 
McCammon  presents  a  formula  for  calculating  die  scattering  strength  based  on  the  surface 
slope  and  rms  slope: 


8l  is  the  local  surface  slope.  6Um  i*  ^  nns  slope  and  is  given  by 


(3.4) 


where  is  die  rms  surfiare  roughness  and  I.  is  the  spatial  onrelation  length.  A  is  a 
scaling  factor  which  varies  depending  on  the  surface  conqiosidon.  This  model  has  been 
successfully  applied  in  the  modeling  of  forward  scatter  and  has  been  included  in  the 


BISSM2  Modified  Bistatic  Scattering  Strengdi  Model.  (McCammon,  1993) 


This  model  qipears  quite  similar  to  the  low-frequency  scattering  function 
for  a  rough  surface  with  an  isotropic  Gaussian  COTrelation  function  developed  by  Tolstoy 
and  Clay  (1987).  This  function  oily  considers  the  incoherent  scattered  radiation  relative 


to  the  incident  radiation: 
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Y»«»  f*(»  (ny 


(3J) 


91  is  the  Rayleigh  reClectioQ  confident  of  Ae  surfooe,  and  for  an^es  less  than  critical. 


has  magninide  equal  to  1.  y  ^  vertical  wave  number  difGnence  for  the  incident  and 
scattered  rays  and  is  given  by 


Y—^ooa6|-^OQaSJ  (3.6) 

where  0i  and  6,  •>'6  die  incident  and  scattered  angles  relative  to  the  vertical  and  I;  is  the 
wave  number  given  by 


(3.7) 

c 

depending  <»  £requracy/in  Hertz  and  speed  of  sound  c  in  meters  per  second,  o  is  the 
ims  rouftiuiess  of  die  surface  in  meters.  f(0)  is  also  an  angle  dependent  tenn  which, 
assuming  a  two  dimenaional  problem,  becomes; 


f(e)« 


1  -^coaSiCoaBj-sindiSiiiSj 
ooa6,'t.ooa02 


(3.8) 


r/  is  approximately  the  correlation  length  r,  of  the  botuxn,  such  dutt 


ri  iP 


(3.9) 


where  R  range  of  the  receiver.  In  the  far  field,  r/  »  Finally  k  is  the  horizontal  wave 


number  difference,  given  in  the  following  equation: 
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(3.10) 


K  -  JbCdaOj-dnfij) 

Substituting  these  definitions  into  the  original  equation,  we  see  tite  somewhat  more 
familiar  focm: 

Sy-lO  logj  (i)  k*6*  [l+coa(e,+ej)f  e'^ 

CmiqMred  to  McCammon’s  fonnula  for  facet  scattering  strengdi,  we  see  an  exptmential 
tenn  whkh  depends  on  the  incident  angle  of  radiation  and  the  spatial  correlation  length 
of  the  bottom.  The  A  parameter  in  McCammon*s  formula  appears  to  be  equal  to  the 
remainder  of  the  equation.  These  two  equations  (33  and  3.11)  are  significant  because 
they  were  to  be  used  to  check  whether  the  backscattering  strength  model  used  in  the 
prediction  routine  is  consistent  with  the  approximate  bottom  roughness  in  the  Barents  Sea. 

b,  Bragg  Scatter 

Bragg  scatter  is  qiplicable  to  low  grazing  angles,  where  ccunposite 
roughness  ctxnputes  scattering  based  on  interacticms  with  resonant  components  in  the 
roughness  of  the  surface.  Hence,  for  those  cmnponents  (facets)  to  be  significant,  they 
must  have  spatial  wavelengths  equal  to  twice  that  of  the  incident  sound.  This  concept  of 
resonances  is  similar  to  that  seen  in  diffiactitxi  gratings,  hence  the  name  Bragg  scatter. 
The  input  parameters  for  this  model  are  the  surface  height  power  spectrum  W(K)  at  the 
resonaiKre  wavenumber  K  ^  k  (Eq.  3.10),  which  can  be  derived  fnmi  the  surface  rms 
roughness  and  correlation  length,  given  the  surface  roughness  has  a  Gaussian  distribution. 
Bragg  backscattering  strength  in  two  dimensions  is  given  by 
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5^.  unJ*" 


(3.12) 


where  the  angles  6«  and  6j  are  die  grazing  angles  oi  the  incident  and  scanned  ener^,  as 
in  Eq.  3.3  (McCammon,  1993).  Again,  dns  dieoty  dqiends  on  knowledge  of  the  assumed 


wavmundier  dependence  die  qiectrum  W.  which  in  general  is  ntM  well  documented  for 
most  of  the  ocean  bottoms  on  the  planet  This  model  is  based  solely  on  a  statistical 
model  oi  the  physical  structure  of  die  reflecting  surface.  Hence,  effects  due  to 
deterministic  structures  (bottom  features),  vdume  scattering  effects  (subbottom  scattering), 
or  refraction  (the  basis  of  discrete  backscaner  models)  are  not  accounted  for.  Asaresult 
the  effective  fourth  power  dependence  on  die  sine  of  the  grazing  angle  causes  to  be 
inaccurate  (low)  below  20”. 


2.  Lambert’s  Law 

The  second  theory  is  Lamben's  law  of  diffuse  scattering.  As  noted  by  Urick, 
Lambert’s  law  is,  "a  ^pe  of  angular  vaiiadtm  which  many  rough  surfaces  appear  to 
satisfy  for  the  scattering  of  both  souml  and  light”  (Urick,  1983)  Noting  the  original 
q^licaticm  of  this  "law”,  he  continues: 

Although  many  materials  follow  Lambert’s  law  closely  in  scattering  light,  none 
does  so  exacdy.  Lambert’s  law  tqiplies  specifically  to  the  radiation  of  light  by 
radiant  absorptive  r.aiterials;  the  "law”  should  properly  be  called  Lambert’s  "rule" 
for  scattering.  Nevertheless,  it  is,  as  we  have  seen,  a  good  description  of  the 
backscattering  of  smind  by  very  rough  bottoms.  (Urick,  1983) 

The  original  triplication  of  Lambert’s  law  was  to  optics,  in  which  it  is  known 
as  Lambert’s  cosine  law  of  diffuse  reflection.  As  described  by  Houstoun  (1938),  a 
portion  dS  of  a  surface  has  apparent  bri^tness  which  varies  with  the  sine  of  the  grazing 
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dS 


Figure  3^:  Geometiy  for  discussion  of  Lambert’s  law  of  diffuse  reflection. 


angle  of  the  incident  light,  but  which  does  not  vary  with  the  grazing  angle  of  the  reflected 
light  Specifically, 


8in(0j) 


(3.13) 


where,  as  shown  in  Hgure  3.4,  pi  is  a  material  dependent  ccmstant,  1.—^-.  is  the  intensity 
oi  the  inddem  light,  6f  is  the  grazing  angle  of  the  incident  light  and  r  is  the  distance 
between  dw  suiftce  and  the  observer.  ^  ^  sin(0|)  is  the  total  light  received  by 
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the  aiea  tff  due  to  light  of  intensity  /.  A  pnctical  exanq>le  ai  diis  relationship  is  a 
luminous  sphere,  for  example,  a  red-hot  copper  sphere  heated  over  a  flame.  When 
viewed  firom  a  distance,  it  i^peais  as  a  disc  of  uniform  tightness  over  the  entire  surface 
when  viewed  from  any  angle.  Hence,  the  radiation  firom  an  element  of  the  surface 
noaking  an  angle  6,  with  that  element  of  ti»  surface  must  ctmsequendy  be  proportional 
to  sin(6«)  (Houstoun,  1938).  The  name  "cosiite  law”  comes  from  the  fact  that  in  its 
original  develqmient,  angle  were  measured  with  the  normal  to  the  surface,  not  the  grazing 
angle  with  the  surface  itself. 

In  the  application  of  Lambert's  law  to  acoustics,  the  scattering  stienvth  varies 
as  the  product  of  the  sine  of  the  incident  aiul  scattered  grazing  angles  times  a  constant  p 
that  is  a  function  of  the  material  pr(^)erties  of  the  scattering  surface.  Thus,  scattering 
strength  is  given  by: 

Sf^ni^^lOlogCpaiiiBgSinfl^)  «  p^-^lOlogCsinOoSiiiO^)  (3.14) 

If  all  of  the  incident  energy  were  redistributed  in  the  water  column,  with  none  lost  by 
transmissicMi  into  the  botumi,  it  can  be  shown  by  integration  that  p  =  l/n,  or  lOlog(p)  = 
p^  s  -4.97  dB.  In  1961,  K.  V.  Mackenzie  reptmed  success  using  Lambert’s  law  with  p^t 
s  .29  dB  to  model  bottom  backscattering  in  deep  water  at  530  and  1030  Hz.  His 
conclusions  are  not  viewed  as  a  rigorous  theoretical  result,  but  instead  a  simple 
geometrical  argument  which  fits  the  angle  dependence  of  ”a  lot"  of  measured  bottom 
backscattering  data  at  moderate  grazing  angles,  specifically  20°  to  60°.  (Mackenzie, 
1961) 
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Figure  3J5:  Reverbendoo  fidlowing  a  2  lb.  e)q)losive  charge  detonating  at  244 
m  in  water  1^0  m  deep,  as  observed  with  a  nearby  hydrc^hone  at  a  depth  of  41 
m.  Filter  band  1  to  2  Idlz.  (Urick,  1983) 

The  choice  of  20"  to  60"  as  limits  fimr  the  implication  is  not  arbitrary.  Instead, 
they  result  fimn  the  geometry  of  the  eiqieximents  that  were  conducted  in  deep  water  and 
the  method  by  which  the  results  were  analyzed.  As  shown  in  Hgure  3.5,  reverberation 
due  to  the  first  surface  reflections,  volume  scattering,  scattering  in  the  deep  scattering 
layer,  and  bottom  backscatter  are  separated  in  time.  At  improximately  2.4  seconds,  the 
first  bottom  returns,  the  so-called  fathometer  returns^  come  in.  At  2.7  seconds  the 
fatlxHneter  return  fiom  the  surface  reflectitm  of  the  explositm  arrives  at  the  receiver, 
followed  by  a  Imig  (~2  sectuid)  bottcmi  reverberation  return.  At  2.7  seconds,  the 
improximate  grazing  angle  of  the  energy  interacting  with  the  bottom  is  60".  At  4.8 
sectmds,  after  v^iich  the  next  set  of  bottom-surface  reflections  arrive  at  the  receiver,  the 
bottmn  grazing  angte  of  the  direct  path  energy  (energy  which  has  traveled  through  the 
water  column  only,  with  no  interactions  with  the  surface  or  bottom)  has  decreased  to 
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^fpnwkniiely  20*.  Hence,  die  limittfbr  die  ^fdkal^ty  of  ihbifionofiiatfcreveiberBtion 
measureniem  arise  naturally  ftom  the  experimental  prooedure.  The  tenn  monostatic  is 
q>proi»iate  in  this  context  because  the  only  scattered  padis  wdiich  are  significant  are  those 
which  follow  the  same  grazing  angle  as  the  incident  energy. 

The  bistatic  dependence  of  Lambert's  law  has  been  examined  by  P.  B. 
Schmidt,  again  in  deep  water,  in  a  1971  piq;)er.  He  found  good  agreement  with  bottom 
scattering  data  for  incident  angles  fiom  3.5*  to  78*  with  scattered  return  angles  from  21* 
to  84*  (McCammon,  1993).  The  value  of  Pus  must  be  set  empirically,  however,  which 
requires  considerable  analysis  in  finding  die  variation  of  with  frequency,  azimuth 
angle,  and  sediment  roughness.  This  theory  has  been  successfully  implemented  by  J.  W. 
Caruthen,  et  aL  in  the  Modified  Bistadc  Scattering  Strength  Model  (BISSM2) 
(McCammon,  1993).  They  found  -3U  dB  to  best  fit  their  data. 

In  shallow  water,  the  problem  becomes  somewhat  mote  complicated.  Hgure 
3.6  shows  the  return  from  a  1.8  lb  char^  used  in  the  BSPFEX.  Instead  of  having  clearly 
separated  volume,  surface,  and  bottom  returns,  they  instead  all  come  in  simultaneously. 
The  implications  of  this  are  two-fold.  First,  the  dominant  scattering  mechanism  must  be 
used  to  predict  the  leverbetatitMi  level  (in  this  case,  bottom  reverberation).  Second,  the 
problem  is  truly  bistadc  in  nature,  since  dw  travel  time  difference  between  direct  path 
rays  and  rays  which  have  multiple  interactions  with  the  surface  and  bottom  are  small,  on 
the  order  cd  the  decay  time  of  the  bubble  pulse  of  the  explosion  (see  Section  C  for  a 
discussion  of  the  bubble  pulse  and  <xher  underwater  explosion  phenomena).  As  a  result, 
for  each  bottom  interaction,  there  are  multiple  scattered  angles  which  significantly 
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IMt  82  SUS  Uvti  at  1  ^  km..  Racaivsd  at  Phon*  0  of  Array 


Figure  RBveri)eratioa  following  a  1.8  lb  explosive  charge  (Mk  82  SUS) 
detonating  at  91.4  m  in  water  275  m  deep,  as  observed  with  a  hydiophtnie  at  a 
depth  of  121  m.  Filter  band  50  to  500  Hz. 

amtribnte  to  the  reverberation.  Only  <xie  of  these  multiple  scattered  angles  is  the 
monostatic  return;  the  scattering  is  dius  bistatk;  in  elevation.  As  will  be  shown  in  Chapter 
rv,  the  scattering  is  in  fact  bistatic  in  azimuth  also. 

In  Ouster  IV,  the  method  used  to  predict  botttnn  reverberation  in  the  Barents 
Sea  will  be  based  on  an  io^lementation  of  Lambert’s  law  in  shallow  water  taking 
advantage  of  tte  vertically  bistatic  nature  of  .  _  This  problem  could  not  be  attacked 
in  diis  manner  in  the  past  because  of  the  lack  of  adequate  cmnputing  power. 
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The  tenenl  fom  of  the  omnidhectioiuJ  scattering  theory,  also  known  as  the 
Lommel-Seeliger  law,  is  (McCammon,  1993): 

5^^*10]og(|iaiiii^)  «  |i^-^lO]Qg(rin00)  (3.1S) 

Again,  ii^  is  dqpendent  ujK»  the  bottom  sediment  properties,  fiequency,  and  azimuth  (of 
the  scattered  energy).  This  theory  differs  from  Lafid)e(t*s  law  in  that  it  only  depends  on 
the  sine  of  die  incident  energy*s  grazing  angle.  McCammon  (1993)  argues  that  at  low 
grazing  angtes  (less  dian  30”),  bodi  monostatic  and  bistadc  data  support  the  first  power 
dependence  on  8in(6^.  This  theory  wUl  be  tested  in  Chapter  IV. 

C.  UNDERWATER  EXPLOSION  THEORY 

Underwater  explosions  have  been  used  for  many  years  as  acoustic  source  in  many 
diverse  applications  (Western,  1960): 

•  Geophysical  jnospecting,  especially  by  the  oil  industry  (where  their  use  has  been 
quite  successful). 

•  Underwater  signaling,  including  distress  signaling. 

•  Sound  sources  for  explosive  echo  ranging  against  submarine  targets  (though  such 
use  is  not  popular  amraig  submariners,  including  the  author). 

•  Underwater  sound  propagation  experiments,  including  reverberation  studies. 

In  die  latter  ai^lication,  they  have  become  the  de  facto  standard  source  due  to  their 
energy  density,  near  impulse-like  characteristics,  and  low  cost  They  also  have  distinct 
advantages  over  odier  impulse-like  sources  (such  as  spaikers)  in  that  they  are  easy  to 
deploy  and  are  not  limited  in  depdi.  This  is  not  to  say  they  are  perfect,  however,  since 
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the  signal  they  (noduoe  is  somewhat  random  and  is  riddled  with  distinguishing 
characteristics  which  vary  dqiending  on  die  type  of  explosive/detonator,  the  fusing 
mechanism,  the  size  of  the  charge,  and  the  deitth  of  detonation.  Since  diey  have  been 
used  for  so  long,  many  of  diese  characteristics  have  been  measured  and  quantized,  and 
hence  diey  have  become  relatively  predictabte  sources  of  sound. 

L  The  Underwater  Explosion  and  Gas  Globe 

When  a  high  eiqilosive  is  detonated,  the  detonation  wave  pix^iagates  through 
the  mass  of  plosive  material,  converting  it  to  an  incandescent  gas  at  a  very  high 
jnessure.  In  hi^  eiqilosives  like  TNT,  this  detonation  wave  or  "shock"  wave  is 
supersonic  with  velod^  between  5,000  and  10,000  m/s.  A  spherically  symmetric  shock 
wave  is  then  radiated  into  die  surrounding  water  with  a  pressure  rise  that  is  practically 
instantaneous,  followed  by  an  exponential  decay  with  a  time  constant  on  the  order  of  a 
fraction  of  a  millisecond.  The  gas  bubble  or  "globe"  created  by  the  explosion  grows  with 
the  shock  wave,  but  doe  to  inertia  it  overshoots  its  equilibrium  radius,  and  hence  the 
radiated  pressure  becomes  slighdy  negative  as  the  bubble  begins  to  contract  As  it 
contracts  it  again  overshoots  its  equilibriiun  radius,  and  when  it  reaches  its  minimum 
volume,  a  positive  pulse  with  strength  and  duration  conqiarable  to  the  initial  shock  wave 
is  develqped.  Successive  pulses  grow  increasingly  weaker.  This  danqied  radial 
osdllatitm  cmitinues,  creating  positive  pulses  at  each  point  of  minimum  radius,  until  the 
globe  either  reaches  the  surface  or  breaks  up  into  smaller  parts.  The  oscillatoiy  signal 
of  these  explosions  is  known  as  the  bubble  pulse.  The  pressure  signature  of  an 
underwater  explosicm  at  short  ranges  consists  of  the  shock  wave  followed  by  a  small 
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muwtMr  dew  Am  10)  of  bidAle  puliec.  Ai  loafer  nuifes  or  in  sliallow  water,  the 
received  lignal  is  oonpUcated  by  miitripath  effects  which  create  more  positive  and 
negative  ioofNilses  for  each  path,  with  suifece  reflections  causing  negative-going  inq>ulses. 
(Weston.  1960.  and  Urick.  1983) 

2.  The  Bubble  Pidae 

The  time  difference  bMveen  the  shock  wave  and  die  first  bubble  pulse  is 
determined  by  the  weight  of  the  charge,  the  type  of  explosive,  and  die  depth  at  which  it 
is  detonated.  This  time  difference  T  is  inedictabte  and  can  be  found  theoretically  (Urick. 
1983): 


Kw^ 

*  (J+33)^ 


(3.16) 


K  depends  on  the  type  of  explosive  (JT  »  4.36  for  TNT),  w  is  die  charge  weight  in 
pounds,  and  dis  the  detonation  depth  in  feet  The  reciprocal  of  T  is  known  as  the  bubble 
pidse  frequency  and  is  measurable  fiom  data.  For  a  Mk  82  SUS  charge  detonating  at  300 
ft.  the  predicted  bubble  pulse  frequency  is  23  Hz.  whidi  agrees  well  with  the  Barents  Sea 
data.  This  measurement  is  inqiortant.  since  it  proves  that  the  SUS  charges  were 
detonating  at  the  correct  depdi  widi  the  conect  strength  in  the  shock  wave. 

Eneigy  is  radiaied  as  long  as  the  bubble  continues  to  oscillate,  but  fOT  large 
charges  and  shallow  depdis  (i.e.  whoi  the  bubble  is  large),  the  bubble  rises  iqipieciably 
during  one  oscillation,  resulting  in  a  new  source  depth  and  a  different  multipath  problenL 
Additionally,  it  affects  the  energy  spectrum  of  the  source,  since  the  frequency  of 
oscillation  will  shift  as  the  babble  migrates  toward  the  surface.  Because  of  the  possible 
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inylkariom  of  diis  btAUt  mignuioiu  it  is  imponant  to  know  when  the  effect  is 
significttL  Weston  (1960)  proposes  that  the  critical  depth  below  which  this  effect  is 
small  is  given  by 

<4  -  (3.17) 

wheie  W  is  the  charge  weight  in  pounds.  This  coneq)onds  to  232  ft  for  a  1.8  lb  charge. 
Since  die  Mk  82  SUS  detonated  at  3(X)  ft,  Inibble  migration  is  not  a  concern  in  the 
BSnPEX. 

3.  The  Encrgy/Fowcr  Spectral  Density 

Since  the  pressure  signature  of  a  SUS  is  by  its  very  nature  random,  it  is 
necessary  to  move  to  the  ftequency  domain  in  order  to  avoid  complications  caused  by 
modeling  such  an  eiqilosion  in  the  time  domain.  This  riiift  also  entails  a  choice  of  either 
energy  flux  densUy  or  power  spectral  density  ai  the  shock  wave  and  first  few  bubble 
pulses  (which  contain  almost  all  of  die  energy  generated  by  the  explosion).  This  thesis 
will  use  power  qiectral  densities  in  units  ci  decibels  referenced  to  1  pPS?/Hz  @  1  meter. 
Weston  and  Urick  both  provide  qiiectra  ftr  charges.  Weston’s  work,  while  excellent,  is 
someudut  dated,  Urkk  indudes  a  diaclaimer  with  his  spectrum  calculadcm  in  which 
describes  his  curves  are,  "based  tqxn,  and  are  useful  only  for,  broad-band  tevels,  and  are 
misleading  for  narrow-band  work  in  not  riiowing  the  fine-frequency  spectral  variation 
caused  by  die  multipulse  character  of  an  explosive  signal."  (Urick,  1971)  Hence,  a  more 
recent  source  is  needed. 
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Figure  3.7:  CoaqMtison  erf*  Mk  82  SUS  source  level  q)ectn. 


Subsequent  lo  die  above  reseaichen,  source  level  measumnents  for  various 
types  of  charges  were  made  by  Stockhausen  (1964).  Christian  and  Blaik  (196S),  Turner 
and  ScriiQger  (1S170).  Gaqrin  and  Shuler  (1971X  tuxi  Chapman  (1987).  Prior  to  Chapman, 
the  moat  widely  used  values  for  the  source  levels  of  shaUow  SUS  were  Gaspin  and 
Shuler's.  Their  estimatBS  were  derived  fipom  hand-drawn  waveforms  that  were  based  on 
experimental  measurements  of  the  waveforms  of  deep  eiqilosions  extrapolated  to  shallow 
dqidis.  Between  Gaspin  and  Shuler's  work  and  Quqmuui's,  there  were  several 
measurements  peifonned  that  were  in  general  agreement  with  the  Gaspin  and  Shuler 


fwedictioiis,  but  there  were  considenble  differences  among  the  dau  sets.  (Chapman. 
1987) 

Chapman’s  results  are  presented  in  1/3  octave  bands  from  5  to  630  Hz.  At 
least  30  charges  (1.8  lb  SUS)  were  avenged  for  each  detonation  depth  reptxied.  The 
standard  deviati<»s  were  within  1  dB  for  mott  frequency  bands  (from  SO  to  SOO  Hz  all 
standard  deviations  were  less  dtan  1  dB).  Chainnan  notes  that  actual  explosim  depths 
were  consistendy  deep  in  all  measurements.  He  theorized,  with  suppcn  frmn  a  previous 
experiment,  that  this  was  due  to  the  relatively  cold  water  (9**  C)  in  which  the  experiment 
was  conducted.  Instead  of  the  charges  averaging  the  nominal  91.44  m  (300  ft),  he  noted 
an  avenge  depth  of  99.6  m  (326  ft).  This  effect  should  be  insignificant  in  the  BSPFEX, 
as  the  effect  bubble  pulse  frequency  is  sli^t  Also,  since  the  propagation  model 
assumes  a  monostatic  problem,  a  deeper  detonation  actually  moves  the  detonation  closer 
to  the  top  of  the  airey,  making  the  mcmostatic  assumption  more  viable. 

Quqnnan’s  results  agree  well  with  Gaspin  and  Shuler’s,  with  the  possible 
exception  of  the  lower  frequencies.  Quqnnan  points  out  cleariy,  however,  that  this  effect 
was  expected  because  Gaspin  and  Shuler  did  not  specifically  account  for  the  impulse  of 
the  waveform,  which  would  throw  off  their  low  frequency  predicticms  (Chapman,  1987). 
In  this  thesis,  the  source  levels  measured  by  Chapman  will  be  used,  utilizing  interpolated 
values  (cubic  spline  fit)  to  obtain  an  approximate  source  level  every  12.5  Hz  from  50  to 
500  Hz.  This  is  shown  by  the  dashed  line  in  Hgure  3.7. 
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IV.  METHOD  AND  RESULTS 


A.  PREDICTION  METHOD 

As  described  in  Qiapta  I,  propagatioa  modeling  for  the  Barents  Sea  has  been 
perfonned  using  lay  traces  provided  by  HARPO.  The  propagation  naodel  from  HARPO 
is  then  used  along  with  die  SUS  source  model  described  in  Chapter  m  and  either 
Lambert's  law  or  oauiidirectional  scattering  bottom  backscatteiing  functions  to  produce 
a  predicted  reverberation  envelope.  Each  oi  the  components  in  the  prediction  routine  are 
described  in  die  sections  below. 

L  HARPO 

HARPO  (Hamiltonian  Acoustic  Ray-tracing  Program  for  the  Ocean)  is  a 
FORTRAN  cooqiuter  program  that  traces  the  three-dimensional  paths  of  acoustic  rays 
through  model  oceans  by  numerically  integrating  Hamilmn’s  equations,  which  are  a 
differential  expresskm  trf  Fermat’s  principle.  R  was  developed  by  R.  Michael  Jones,  J. 
P.  Ril^,  and  T.  M.  Georges  at  the  Wave  Prt^gation  Laboratory  in  Boulder,  Colorado 
(an  Environmental  Research  Laborainy  of  the  National  (Jceanic  and  Atmospheric 
Adnnninration).  As  described  in  Jones,  et  aL  (1986),  it  uses  continuous  models  for  the 
stnind-speed  distribution  and  bathymetry  in  die  ocean,  hence  avoiding  the  false  caustics 
and  diacontinuous  raypath  (Rtqierties  encountered  in  conventional  ray-tracing  methods. 
Only  gecxnetrical  effects  are  accounted  for  in  the  ray-trace;  diat  is,  no  diffraction  or 
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IMvdat  reflection  GOtrectiolu  are  apfrfied  by  HARPO 1^^  Anqilitude  akxig  a  nypuh  is 
not  ejqplidtly  coooputed;  that  is  left  to  a  post-processing  program  which  accounts  for 
boundary  interactions  and  losses  due  to  ray  tube  spreading  and  absorption.  The  program 
does  produce  a  stq>-by-stq)  account  of  each  ray*s  {nogress,  including  all  turning  points 
and  boundary  interactions,  and  outputs  die  data  in  machine-readable  form  for  post¬ 
processing.  Tire  post-processing  tool  chosen  for  this  thesis  was  MATLAB. 

The  inputs  for  HARPO  in  its  present  form  (on  the  Naval  Postgraduate 
School’s  Cray  Y-MP  EL98)  are  a  data  file  which  contains  the  three-dimensional  sound 
speed  field  in  the  region  (sound  speed  at  each  longitude/latitude/depth  point),  a  data  file 
which  contains  the  two-dimensional  bathymetry  for  the  same  region  (depth  at  each 
longitude/latitude  pair),  the  soutce/receiver  location,  die  desired  launch  elevation  angles, 
and  the  desired  launch  azimuth.  It  then  computes  closed-form  expressions  for  the 
bathymetry  and  the  sound  speed  field  using  a  user-specified  number  of  ’’empirical 
orthogonal  functions”  (EOF’s)  to  model  the  sound  speed  field  (Chiu,  et  al.,  1993).  The 
number  of  EOF’s  required  depends  on  the  miginal  number  of  sound  velocity  profiles 
(SVP’s)  entered  and  the  maximum  eim’  desired  in  the  continuous  representation  of  the 
field. 

It  should  be  noted  that  HARPO  cannot  directly  compute  the  raypaths  that 
connect  a  specified  source  and  receiver  (or  scatterer  in  the  case  of  reverberation).  To  find 
such  eigenraySy  one  usually  launches  a  fan  of  rays  at  small  increments  in  elevation  and/or 
azinmth  angle  and  linearly  interpolates  for  the  rays  that  reach  the  desired  receiver 
locatitm.  In  the  specific  case  of  bottmn  reverberation,  since  each  bottom  interaction  is 
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Original  (*)  and  Extrapolated  (o)  CTO  Locations  (w/isobaths) 


Figure  4.1:  CTD  locations  and  extnqx)lated  points  along  isbaths. 


identified  in  HARPO’s  ouqtut  "rayset,"  this  turns  out  to  be  an  interpolation  along  the 
bottcmi  only  (assuming  the  bottom  is  locally  flat).. 

The  bathymetry  for  the  Barents  Sea  was  obtained  using  tte  WORLDBATHY 
database  in  the  Naval  Postgraduate  School’s  Oceanogr^>hy  laboratory.  The  data 
originally  came  from  DBDB5,  a  standard  U.  S.  Navy  worid  bathymetric  database.  The 
sound  speed  field  was  measured  as  part  of  the  BSPFEX.  A  total  of  72 
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Conduciiviiy/Tcinperaturei/Dqwh  (CTD)  measurements  were  taken,  56  of  which  were 
within  38  km  of  the  array  and  hence  were  usable  in  the  reverberation  experiment  To 
cover  the  entire  area,  the  CTD  measurements  were  copied  to  new  locations  along  isobaths, 
smoothed  using  an  1 1  point  Hamming  smoother,  and  then  an  inverse  distance  algorithm 
(MATLAB,  1992)  was  q)plied  to  interpolate  the  sound  velocity  field  at  regular  grid 
locations  fw  10  longitudes  and  9  latitudes  around  the  array.  For  the  three-dimensional 
sound  speed  field  in  the  Barents  Sea  expoiment  area,  a  total  of  90  SVP’s  were  used  to 
estimate  the  sound  speed  field.  To  keep  cmx  less  tiuui  0.1  m/s,  26  EOF’s  were  required. 

To  estimate  the  horizontally  bistatic  nature  of  the  experiment  (the  source  SUS 
and  receiving  array  were  ~1  km  rq)art)  using  a  monostatic  propagation  model,  a 
soutceAeceiver  position  tiuit  was  in  the  center  of  the  SUS  drop  points  and  array  location 
was  used.  From  that  point,  rays  were  launched  on  several  azimutiis  over  a  fan  of 
elevation  angles.  The  runs  were  originally  performed  using  a  three-dimensional 
environment,  but  it  was  quickly  discovered  that  doing  so  violated  the  assumptions  made 
in  modeling  the  reverberation  using  an  N  x  2D  approximation.  To  correct  this 
discrepancy,  the  final  HARPO  runs  were  performed  using  two-dimensional  "slices"  from 
the  three-dimensional  environmental  data.  Four  launch  azimuths  (000, 090, 180  and  270) 
were  used,  with  rays  traced  every  0.5®  in  elevation  from  -60®  to  -20®  and  from  -•■20®  to 
-h60®,  and  every  0.1®  in  elevaticm  from  -20®  to  -i-20®.  The  results  of  both  the  2D  and  3D 
data  runs  are  described  later  in  tins  chaprer. 
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2.  Poit-processing  HARPO:  rbreakbr.ni 

rbreakbrm  is  a  MATLAB  script  which  reads  in  the  output  data  rile  from 
HARPO  and  processes  it  to  find  eigenrays  between  the  sourceAeceiver  location  and 
desired  range  points  along  the  azimuth  cm  which  HARPO  was  run.  For  each  eigenray  to 
a  given  botrom  point,  it  determines  the  ray's  amplitude,  path  length,  horizontal  range, 
depth,  phase,  arrival  angle  and  travel  time.  The  local  sound  speed  at  the  bottom,  original 
launch  angle  of  the  ray,  and  vertical  distance  between  it  and  the  adjacent  ray  used  for  ray 
tube  spreading  calculation  are  also  stored. 

Ray  path  lengths  are  determined  by  simply  summing  the  straight  line  distances 
between  HARPO  rayset  ouQtut  points.  Stepping  tiirough  all  ouq)ut  points  for  each  ray, 
boundary  interactions  and  turning  points  are  accounted  for  individually.  At  each,  the  type 
of  interaction  is  identified  and  numbered  using  a  system  described  by  Ranchi  et  al. 
(1984).  Each  interaction  or  turning  point  is  numbered  as  an  order  contour^  with  all 
surface  inteiactions/upper  turning  points  receiving  odd  numbers,  and  bottom 
inteiactions/lower  turning  points  even.  If  the  first  interaction  is  a  bottom  interaction,  it 
is  hence  identified  as  #2.  Each  is  also  ictentified  as  to  the  type  of  point  (boundary 
interaction  or  turning  point).  This  numbering  method  turns  out  to  be  key  in  later 
processing  when  interpolation  is  performed  between  adjacent  launch  angle  rays  at  bottom 
interactions  with  the  same  order  contour. 

Amplitude  calculations  for  a  ray  at  a  given  bottom  interaction  are  determined 
by  the  combined  effects  of  surface  and  bottom  interactions  the  ray  has  experienced  up  to 
that  point  and  ray  tube  spreading.  At  each  surface  interaction,  the  loss  mechanism 
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WiWHiwi  it  idtoi  ID  iCiWBitBg  Ion  only  osing  te  Edctit  ooheient  reflectioii  coefBcient  Rc 
such  dm 


(4.1) 


where  SI,  die  Rayleigh  reflecdoo  ooefBciem,  is  equal  to  1,  is  the  tins  waveheight 
in  meters  (less  than  0.1  meter  during  the  etqperiement;  0  meters  used  in  analyses),  and  y 
is  the  horizontal  wave  number  CTolsioy,  et  aL,  1987).  For  each  bottom  interactitHi,  the 
loss  mwrhaiwMn  asstimed  is  due  to  scattering  loss  using  the  Eckan  coherent  reflection 
coeffidem  along  with  91.  For  ang^  less  than  critical  (68**  for  die  bottom  in  the  vicinity 
of  the  aiTi^.  as  noeasured  using  the  ground  wave  oi  the  SUS  charges)  the  magnitude  of 
the  Rayleigh  lefiectioo  confident  is  1.  Rays  which  have  angles  higher  than  this  are  so 
severely  attenuated  that  dieir  contribution  to  the  reverberation  level  is  insignificant. 
Hence,  the  magnitude  of  the  specularly  reflected  ray  after  a  bottom  reflection  is  also 
detennined  by  Equation  4.1,  where  —  is  the  ims  roughness  of  the  bottom  in  meters 
(Tolstoy,  et  aL,  1987).  -  was  determined  fiom  fathometer  recordings  to  be 

approximately  0.5  meten.  At  die  fiequendes  of  interest  in  the  BSPFEX,  loss  due  to 
absmption  in  the  bottcxn  is  negligible.  Ray  tube  spreading  is  calculated  using 
conservatkn  of  power  widiin  the  ray  tube  given  by  the  foUowing  equation: 


SO 


(4^) 

F  is  the  acoustic  pressure  at  horizontal  lange  r  firom  die  source,  1%  is  dre  acoustic  pressure 
1  meter  from  die  source,  is  the  angular  difference,  in  radians,  between  adjacent 
launch  angle  rays,  A  is  the  vertical  distance  between  die  two  rays  at  range  r,  p  and  p«  are 
the  water  densities  at  die  bottom  point  and  at  the  source  (assumed  equal),  e  and  are  the 
sound  qieeds  at  die  bottom  point  and  at  the  smirce,  0«  is  the  launch  elevation  and  6  is  the 
elevatitm  angle  the  ray  at  die  bottom  interactitm  (Clay,  et  al.,  1977). 


IMhodol  imtfpolaiion  Bsiwmo  Adfsosni  Rcys 


Rang*<lun) 

Figure  42:  Notional  view  of  interpolation  between  adjacent  launch  angle  rays. 
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Once  all  bonoea  imarecricni  have  been  idenrifled  and  eifenray  data  for  each 


has  been  stoied,  the  program  dien  steps  out  in  range  to  the  points  specified  (at  regular 
intervals)  and  finds  adjacent  launch  angle  rays  which  have  bottom  interactions  with  the 
same  order  contour  which  "straddle"  die  desired  bottom  pcnnc  This  is  illustrated 
grqihfoaUy  in  Figure  4J2^  adiere  the  two  scdid  lines  represent  adjacent  launch  angle  rays 
(which  have  bottom  interactions  with  the  same  order  contour)  which  straddle  the  location 
of  a  desired  ouqiut  point  on  the  boQom.  A  linear  interpolation  is  then  performed  in 
range,  assuming  the  bottom  to  be  locally  flat,  to  find  ray  an^litutte,  path  length,  depth, 
phase,  travel  time,  arrival  angle  and  local  smind  speed  at  the  desired  point 

3.  Combining  Eigatrays,  SUS  Model  and  Scattering  Function:  rlmaker.m 
The  goal  of  the  rtmakerjn  rmitine  is  to  estimate  the  reverberation  level 
prediction  algorithm  suggested  by  Franchi,  et  al.  (1984)  and  more  recendy  by  Bucker,  et 
al.  (1993).  They  both  suggest  a  double  integral  over  the  reverberating  area  with  multiple 
eigenrays  cmmecting  each  "patch"  <£4  on  die  bottom  with  the  sourceAeceiver.  In  order 
to  model  reverberation  in  this  manner  using  a  true  three-dimensional  model,  ray  traces 
would  have  to  be  done  from  each  SUS  loooion  over  a  fan  of  launch  elevations  for  each 
azimuth  in  a  360**  fan  surrounding  the  SUS.  A  similar  run  would  have  to  be  done  for 
each  of  the  other  ditee  SUS  and  fw  dre  vertical  line  array.  Then,  rays  connecting  the 
SUS  and  the  vertical  line  array  with  a  given  patch  on  the  bottom  would  have  to  be  picked 
out  of  the  set,  and  the  appropriate  amplitude/travel  times  assigned  to  each.  This  is  a 
fOTmidable  task,  yet  Bucker,  et  al.,  claim  vo  have  successfully  modeled  reverberation  in 
this  way  with  some  sinqilificatimis  in  their  RUMBLE  product  (Bucker,  et  al.,  1993). 
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Confining  the  problem  to  a  twonlimensional  slice  of  the  ocean  and  assuming 
a  horizontally  monostadc  problem  lessens  the  cmiqmtatimial  load  significantly,  and 
enables  the  routine  to  nm  within  a  reasonable  amount  of  time.  It  also  ensures  that  the 
rays  connecting  the  source  to  a  given  scattering  patch  are  the  same  as  die  rays  which 
omnect  that  patch  back  to  the  receiver.  The  reverberation  level  at  die  receiver  due  to  a 
source  with  intensity  l(t)  is  thus 

m-U  E  E  Vf-Tfn-m  " 

The  two  T  terms  represent  the  travel  times  along  the  ith  raypath  from  the  source  to  the 
bottom  point  and  for  the  jth  return  path  in  the  set  of  raypaths  represented  by  the  set  m(r) 
connecting  the  source/receiver  and  the  bottran  patch  dA  which  is  at  range  r.  The  two  B 
terms  reinesent  the  beamforming  factors  to  be  taken  into  account  on  the  source  and 
receivers  and  the  two  L  terms  account  for  losses  in  anqilitude.  a  is  the  scattering  strength 
of  the  bottom,  depending  tm  the  incident  energy  grazing  angle  6/  and  the  scattered 
energy’s  launch  angle  0/.  No  azimudial  dependence  is  modeled.  This  method  is 
displayed  graphically  in  Figure  4.3. 

Ranchi  et  al.  note  that  it  was  not  possible  for  them  to  evaluate  Eq.  4.3 
direcdy  due  to  a  lack  of  a  three-dimensional  prc^gadon  model.  Bucker  et  al.  were  able 
to  evaluate  the  intergral  because  they  had  both  the  model  and  the  required  computing 
assets.  At  the  Naval  Postgraduate  School,  there  were  adequate  computing  facilities  and 
a  three-dimensional  prcqiagadon  model,  but  the  run  time  required  to  properly  model  the 
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Fifure  43:  Olusintiaii  of  leveibentioii  problem. 


reverberatkm  in  shallow  water  proved  prohibitive  at  the  time  of  this  writing. 

Given  die  asssunqitions  stated  above,  and  estimating  die  area  integral  by 
summing  the  contributions  to  reverberation  by  rings  on  die  bottom.  Eq.  4.3  sinqilifies  to 

B  is  now  the  effective  beamwidth  dt  the  point  source/ireoeiver  given  the  regitm  for  which 
the  ray  trace  is  applicable.  For  the  four  azimuth  apprmdmation  to  be  made,  B  »  0.25. 
l(t^  is  die  intensi^  of  the  SUS  at  fiequency/and  is  assumed  ctmstant  at  its  peak  level 
for  0.1  sec.  The  summation  over  k  performs  a  power  sum  of  the  contributions  from  each 
ring  out  to  a  range  dt  ~38km.  placing  the  reverberant  energy  into  the  appnqniate  0.1 
second  time  bin  depending  on  die  two>way,  vertically  bistatic  travel  time  of  die  (/raypath. 
The  two  loss  terms  and  die  bottom  backscattering  strength  are  frequency  dependent 


and  Ti,  lepment  the  imier  and  ouier  radii  of  the  kth  ring  which  has  its  center  at 
and  such  that  rinp  do  not  overly).  Hence,  Ri/tJ)  is  the  predicted 
reverberuiao  peak  envelope  at  the  desired  fiequency. 


B.  BARENTS  SEA  DATA  ANALYSIS 

To  most  doaely  estimaie  a  mnnostatk  problem,  only  die  phone  closest  to  the 
daronation  dqpdi,  phone  0,  the  top  phone,  was  analyzed  for  reveiberation  level. 
Broadband  peak  level  measurements  were  peifonned  using  a  peak  picking  algorithm  and 
a  straight  oonvenkm  fiom  voltage  to  sound  pressure  level  in  decibels: 

ItLCO-  20  kigjivottagfi  X 10^) 

Voltage  levels  were  corrected  for  phone  srosidvity  problems  and  non-zero  mean  as 
described  in  Chapter  n. 

Nanowband  measurements  were  petfoimed  using  a  Short-Time  Fourier  Transfcxm 
(STFT)  power  qiectral  densi^  estimaie  using  ui  AT  >  128  point  rectangular  window,  with 
1/2  window  overhqi  in  segment  intervak  (llienien,  1992).  If  the  fast  Fourier  transform 
(fit)  of  data  vector  x(t)  is  given  by  2r(/),  then  the  narrowband  level  of  x(t)  in  a  1  Hz 
bandwidth  is  given  by 


(4.6) 


is  the  sanqiUng  frequency  (1600  Hz).  The  factm  of  2  is  required  to  maintain 


conservation  of  power  when  only  considering  positive  frequency.  Prediction  results  and 
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C  PREmcnON  RESULTS 

Nanowtand  predictions  were  made  far  1  Hz  bandwidths  oeniered  at  50, 100,  200 

aiid400Hz.  In  the  process  of  developing  the  reauto,  die  hadrscatwaing  coefiBcient  of  the 

Barems  Sea  has  been  deiennine±  Bodi  Landien's  law  and  omnidirectional  scattering 

were  considered,  bi  general,  the  choice  of  the  backscattering  function  tends  to  change 

the  slope  of  the  estimate  with  time,  while  the  choice  of  the  backscanering  strength  shifts 

the  retinaaie  up  and  down  without  slope  change.  Lambert’s  law  provides  a  better  fit  to 

the  measiaed  data.  The  values  of  |i  in  T^le  4.1  were  found  to  best  fit  the  Barents  Sea 

data.  They  are  consistent  widi  the  findings  of  McCammon  (McCammcMi,  1993). 

TABLE  4.1:  BACKSCATTERING  STRENGTH  PARAMEI1ER  pi  FOR 
SEVERAL  FREQUENCIES  IN  THE  BARENTS  SEA. 


Fnqamcj  (Hz) 

50 

100 

200 

400 

pi  (dB  rcmi*) 

>40±15 

-37±10 

-2S±7 

-20±7 

Detennination  oi  which  scattering  functicm  is  better  is  somewhat  subjective,  in  that 
the  chmce  of  parameters  in  the  piqpagation  model  detennine  which  scattering  function 
fits  better.  Hgure  4.4  demonstrates  this  by  showing  the  difference  obtained  using 
Lambert’s  law  and  omnidirectiooal  scattering.  Any  loss  term  which  increases  with  range 
also  changes  the  slope  of  the  estimate,  nxxe  dramatically,  in  fact,  than  the  choice  of  the 
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scattering  function.  Spedfkally.  volume  absoq)tk>n,  suiftce  scattering  loss,  bottom 
scattering  loss,  and  bottom  absorption  are  aU  loss  terms  which  tend  to  act  in  a  roughly 
linear  nsanner  to  increase  transmission  loss  (11.  in  decibels)  with  range.  Therefore, 
perfect,  or  near  perfect,  prc^mgadoo  naodeling  is  required  in  order  to  determine  which 
scattering  function  is  more  accurate.  In  shaltow  water,  this  may  not  be  a  reascmable 
request,  especially  givra  the  nature  of  the  nailtipath  problem. 


Comparison  of  Lambsffs  Law  and  Omnidiraellonal  Seatlaring 


Figure  4,4:  CbtiqMuison  of  Lambert’s  law  with  omnidirectional  scattering  at  200 
Hz,  overlaid  cm  200  Hz  data  from  die  BSPFEX. 


Hgures  4.5  through  4.8  show  die  predictions  and  data  collected  during  the  BSPFEX 
that  generated  the  numbers  in  Table  1.  At  SO  and  100  Hz,  the  propagation  model  is 
suspect  for  two  reastms.  The  rirst  is  die  lack  of  an  adequate  loss  mechanism  to  fully 


•Kplrii  te  dopt  0#  dit  BffFEX  <ba.  TWt  is  eipecialty  \isiUe  at  90  Hz,  is  the 
predictioa  begins  to  increase  at  die  20  second  mark.  The  second  reason  is  the  more 
ftindamwital  question  of  wfaedier  ray  theory  is  conect  at  lower  firequendes.  For  these 
reasons,  die  predictions  at  200  and  400  Hz  are  thou^tt  to  be  much  more  accurate  than 
those  at  the  lower  frequencies.  For  die  validity  of  these  models  m  truly  be  determined 
in  a  general  sense,  a  closer  examination  is  required  using  data  from  several  different 
regions. 

The  trend  indicated  in  die  rising  |t  widi  frequency  is  significant  It  is  also 
consistent  with  rough  surfrice  scattering  theory  as  described  in  Chiyicer  IIL  Referring 
back  to  Equation  3.5,  it  is  clear  that  the  scattering  function  varies  with/in  a  complicated 
fashkm.  From  the  predicticMis  in  Hgures  4.5  through  4.8  and  the  values  listed  in  Table 
4.1,  there  appears  to  be  between  3  and  8  decibels  per  octave  difference  in  scattering 
strength  over  the  frequency  band  considered.  Again,  this  conclusion  should  be  tempered 
with  die  knowledge  that  the  propagatitm  modeling  may  not  be  exact  below  200  Hz. 


25  30 

Tinm(Me) 

Figure  4^:  Reverberation  level  prediction  and  Barents  Sea  results  at  100  Hz. 


Top  view  of  raypalM  launched  on  azimuth  0 


vertical  array  azimuth  OOOT. 


A  significant  effea  not  taken  into  account  in  the  two-^limensional  prediction 
routiiws  is  that  of  horiztmtai  lefiraction.  The  significance  of  this  effect  is  cleariy 
illustrated  in  Hguie  4.9.  In  it  are  shown  the  raypaths  for  three  different  launch  angle  rays 
(+25**,  -140^,  and  4<50*’).  The  bathymetry  proves  to  be  the  dominant  effect  cm  the  rays  as 
they  reftaa  eastward  away  from  the  shallow  region  to  the  nmthwest  of  the  array.  This 
causes  a  two-fold  problem  in  the  reverberadmi  modeL  Rrst,  the  assumption  that  a  ray 
can  be  scattered  to  a  return  path  with  smne  angle  other  than  its  incident  angle  is  invalid, 
unless  a  fan  launch  azimuths  is  used  altmg  with  a  three-dimensional  propagation 


61 


J: 


model.  For  example,  assuimng  that  a  single  launch  azimuth  will  suffice  for,  say.  a  90* 
sector  is  inctmect,  since  this  would  inq)ly  duu  energy  fixxn  a  25“  'aunch  elevation  striking 
the  bottom  at  (74.65*  N,  023.66*  E)  would  smnehow  be  scattered  into  a  45*  received 
elevation  eigenray  which  emanated  from  (74.60*  N.  024.00*  E).  This  is  simply  not  a 
good  assumption.  Of  course,  nxxleling  propagadoi  in  a  three-dimensional  environment 
with  a  two-dinrensional  irxxiel  is  not  a  good  assumption  either,  but  at  least  doing  so  is 
consistent  with  the  method  of  reverberation  jnediction.  The  second  problem  is  the  effect 
of  horizontal  refraction  on  signal  amplitude  due  to  ray  tube  spreading  in  the  horizontal. 
From  Figure  4.9,  this  problem  clearly  appears  to  be  significant.  Hence,  to  accurately 
noodel  the  reverberation  in  a  shallow  regicm,  a  three-dimensional  noodel  should  be  used. 

D.  STATISTICAL  ANALYSIS 

A  statistical  analysis  of  the  BSPFEX  data  has  been  performed.  Since  the 
reverterant  return  from  an  explosion  in  shallow  water  is  by  its  very  nature  a  non- 
stationary  process  in  time,  the  wideband  returns  (sarrqrled  at  16(X)  Hz)  were  broken  into 
5(X)0  point  segments  and  analyzed  individually,  assuming  approximate  stationarity  over 
the  segment.  The  standard  deviation  of  each  segment  was  determined  using 
MATLAB,  and  the  data  segment  was  then  divided  by  o«  to  normalize  the  segment  to 
standard  deviations  of  pressure,  a  unitless  quantity.  This  allowed  segments  to  be 
compared  directly  to  determine  the  underlying  probability  density  function  of  the  entire 
reverberation  signal.  Dividing  the  pressure  by  the  standard  deviation  is  essentially  the 
same  as  dividing  by  the  rms  pressure,  or  normalizing  the  signal  power  to  one  unit  of 
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Figure  4.10:  Nonnalized  probability  density  function  dT  reverberation  signal. 

Dots  show  mean  experimental  values:  vertical  lines  show  95%  confidence 
intervals.  Gaussian  density  overlaid  for  r^erence. 

power.  The  ptobabili^  density  function of  each  segment  was  then  estimated  using 
the  histogram  subroutine  in  MATLAB  (1992)  and  the  equation: 


(4.7) 


where  Jf  is  the  number  of  data  points  used  for  the  estimate.  Ax  is  the  size  of  the  analysis 
bin  in  the  histogram,  and  is  the  number  of  data  points  which  fell  in  the  bin  centered 
on  jt  (Therrien,  1992).  The  average  over  time  was  also  complemented  with  an  ensemble 
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aveiBfe  over  the  four  SUS  returns  to  yield  die  mean  probability  density  function  shown 
in  Figure  4.10.  The  dots  on  the  plot  refnesent  the  mean  probability  density  function 
estimate  over  all  four  SUS  returns,  and  the  vertical  line  through  each  dot  shows  the  95% 
confidence  interval  determined  using 


- 

JC- 


^  I*. 


<  X*  ' 


(4.8) 


where  Jt  is  die  sample  mean  of  the  random  variable  x,  is  the  sample  standard  deviation 
of  X,  is  the  value  of  the  Student  t  distribution  fotn»N~l  (Nis  the  number  of  data 
points  in  the  sample  x),  and  a  =  J  •  P  foraf  P*100  )  percent  confidence  interval. 
(Bendat,  et  al.,  1971) 

Qearly,  the  normalized  wideband  SUS  returns  have  a  stationary  Gaussian 
distribution  (i.e.  a  normal  density  functitm).  Individual  unnmrmalized  segments  are  still 
Gaussian  processes,  but  ate  not  stationary  since  they  have  variance  proportional  to  the 
instantaneous  power  of  die  return  signal  The  erratic  estimate  in  the  vicinity  of  zero 
standard  deviatitms  is  due  to  quantization  noise  of  the  A/D  converter  in  the  measurement 
system.  This  noise  becomes  significant  in  die  normalized  data  late  in  a  given  SUS  return 
when  the  data  is  just  barely  above  background  noise. 

In  order  to  examine  the  statistics  of  narrowband  data,  the  data  were  again  segmented 
in  order  to  cturqiute  fast  Fourier  transfmms  (ffi*s)  of  each  128  point  segment.  In  order 
to  prove  that  each  frequency  bin  has  a  bivariate  Gaussian  distribution,  the  real  and 
imaginary  portions  must  be  shown  to  both  have  Gaussian  distributions.  This  was  done 
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Nonnalizwi  PrababiMy  Oansily  Function  EstimatM  for  Narrowband  Data 


RaaHfftIRL(t)S.  SO  Hz  lmag{fft(RL(t)D.  SO  Hz 


•4  -2  024  ■4  -2  024 


X  X 

Figure  4.11:  Nairowband  probability  density  function  estimates  at  SO,  100,  200, 
and  400  Hz.  Dots  show  mean  experimental  values;  vertical  lines  are  95% 
confidence  intervals.  Gaussian  density  function  overlaid  for  reference. 
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ttsilig  a  mttliod  rimilar  to  thtt  done  on  the  wideband  data  described  above.  The  results 
of  this  analysis  are  shown  in  Hgure  4.11.  As  before,  the  dots  rqireaent  the  sample  mean 
of  the  probability  drasity  function  estimate  over  all  four  SUS  charge  returns.  The  vertical 
lines  are  the  9S%  confidence  intervals.  The  narrowband  data  again  clearly  has  a  Gaussian 
distribution  at  all  finequencies. 

These  results  are  significant,  because  verifying  the  signals  to  be  Gaussian  random 
processes  allows  a  great  deal  of  existing  perfonnance  predictkm  algorithms  to  be  iqiplied 
in  the  analysis  of  active  sonar  systems  in  the  Barents  Sea.  The  results  are  not  surprising, 
since  the  data  cmnes  from  a  multipath  acoustic  problem  which  is  by  nature  randmn  due 
to  the  effects  of  the  oceanogtrq)hic  and  bathymetric  features  over  the  entire  azimuth. 
Nonetheless,  they  are  quite  valuable  to  a  potential  researcher  at  active  sonar  system  model 
for  the  Barents  Sea.  This,  however,  does  not  necessarily  imply  that  these  results  are  true 
in  general.  In  order  to  extend  this  assumption  to  other  shallow  regions,  a  similar  study 
must  be  done  for  that  region. 


V.  CONCLUSIONS 


The  objectives  ci  diis  thesis  were  to  consider  several  methods  of  reverberation 
prediction  and  compare  these  mediods  to  die  results  from  the  BSPFEX,  to  determine  the 
reverberant  characteristics  of  the  Barents  Sea,  and  determine  the  statistical  characteristics 
of  the  reverberation  measured  in  die  BSFFEX.  The  reverberation  predicdon  routines 
which  were  developed  for  diis  shallow  environment  udliaed  muldpath  propagation 
modeling  akmg  with  a  vertically  bistatic  scattering  function  for  the  bottom.  Titty  have 
been  used  to  accurately  model  die  reverbetaiion  due  to  SUS  charges  drcqiped  during  the 
BSPFEX.  It  has  been  determined  that  Landiert’s  cosine  law  accurately  describes  the 
bottom  backscattering  function  using  a  parameter  which  varies  depending  on  bottom 
composition  and  roughness.  Fiirthennore,  diat  parameter  has  been  determined  at  several 
frequencies  for  the  Barents  Sea.  Omnidirectumal  scattering  theory,  howevn,  has  not  been 
ruled  out  as  a  viable  reverbenuion  predictitm  tool.  Finally,  the  statistics  underiying  the 
revoberant  returns  in  the  shallow  water  of  the  Barents  Sea  have  been  proven  to  be 
Gaussian  disoibudons,  both  for  broadband  and  narrowband  data. 

These  results  do  not,  by  themselves,  provufe  all  the  information  required  to  predict 
the  performance  of  a  LFAA  system  in  shallow  water.  They  do  however,  provide  key 
infOTmation  which  is  necessary  to  achieve  dus  goaL  Knowledge  of  which  scattering 
functimi  to  use  and  the  parameters  which  drive  the  scattering  is  essential  to  predict 
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H^iMteriewMtentioiiisaliiiiii^  tf  U  is,  kaov^edfe  of  the  statudcai  propenies 

of  the  reverbenut  retiiint  is  rssmrisi  in  such  assessments. 

Further  research  effon  should  be  dhected  at  modeling  the  reverboation  in  different 
shaltow  regions  using  different  types  of  sources  and  receivers.  Such  modeling  should 
consider  three-dimensknal  effects  bodi  on  piopagatioa  and  scattering. 
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